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Abstract 

We study two-legged spin-1 ladder systems with D 2 xcr symmetry group, where D 2 
is discrete spin rotational symmetry and a means interchain reflection symmetry. The 
system has one trivial phase and seven nontrivial symmetry protected topological (SPT) 
phases. We construct Hamiltonians to realize all of these SPT phases and study the 
phase transitions between them. Our numerical results indicate that there is no direct 
continuous transition between any two SPT phases we studied. We interpret our results 
via topological nonlinear sigma model effective field theory, and further conjecture 
that generally there is no direct continuous transition between two SPT phases in one 
dimension if the symmetry group is discrete at all length scales. 
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1. Introduction 

Symmetry protected topological (SPT) phases are formed by gapped short-range- 
entangled quantum states that do not break any symmetry [1]. Contrary to trivial sym¬ 
metric states, a nontrivial SPT state cannot be transformed into direct product state (or 
Slater determinant state for fermions) of local atomic basis via symmetric local unitary 
transformations. Nontrivial SPT states are characterized by their gapless or degenerate 
edge states at open boundaries, which are protected if the symmetry is reserved. On 
the other hand, when the symmetry is explicitly broken by perturbations, the nontriv¬ 
ial SPT phases can be adiabatically connected to the trivial phase and the boundary 
can be gapped out (for example, see [3]). The Haldane phase [2] in spin-1 chain is a 
one-dimensional (ID) bosonic SPT phase protected by D 2 spin rotation symmetry, or 
time-reversal symmetry, or spatial inversion symmetry [3, 4]. ID bosonic SPT phases 
with onsite symmetry are classified by the projective representations of the symmetry 
group [5, 6, 7]. In higher dimensions, bosonic SPT phases are classified by group co¬ 
homology theory [8], which is consistent with the continuous nonlinear sigma models 
[9] or Chern-Simons theory (in two dimensions)[10]. According to the classification 
theory, new ID SPT phases other than the Haldane phase have been constructed in spin 
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Figure 1: Schematic phase diagram for the Hamiltonian (6). There are four SPT phases, one of which is 
trivial. The SPT phases a, b, c stand for t xz , t y and to phases respectively. The symmetry breaking phases 
Mi, M 2 , Ni, N 2 have order parameters |{5 !,+ )|, |{5 !,_ )|, |(S' 3- )| and |}| respectively. For details 
see Sec. 3. 


chain/ladder models [11, 12, 13]. It was also shown that different SPT phases can be 
distinguished by different responses of their edge states to external symmetry breaking 
fields. In 2D, the boundary of a nontrivial SPT phase is either gapless or symmetry 
breaking [14, 15]. Nontrivial 2D SPT phases, such as integer bose/spin quantum Hall 
phases[16, 17, 18, 19, 20, 21, 22] and bosonic topological insulators [23, 24], have 
been realized. Important progress has also been made in classification and realization 
of SPT phases in three- or higher dimensions [25, 26, 27, 28, 29]. 

In the present work, we will study the model realization of ID SPT phases in spin-1 
ladder systems and the phase transitions between these phases. The symmetry group 
we are considering is an Abelian discrete group D -2 x a. Here the D 2 subgroup is 
a discrete spin rotation symmetry D 2 = {E, R x , Ry, R z } and <7 = {E,P} is the 
interchain exchange symmetry, where E is the identity operation, R x (R y . R z ) is a 
spin rotation of angle 7r along x(y. z) direction, and P is interchain reflection. Since 
T-L 2 (D 2 x a, E/(l)) = Z|, there are eight classes of projective representations, one is 
trivial and the others are nontrivial. Consequently there are seven distinct nontrivial 
SPT phases. The three root SPT phases according to the Ti\ classification are called 
t x , ty, t z , respectively, and all the other phases t xy = t x x f 9 , t xz = t x x t z , t yz = 
t y x t z , to = t x x t y x t Z) I = t x x t x = ty x ty = t z x. t z can be generated 
by stacking two or more of the root phases (where the to phase is the usual Haldane 
phase and I is the trivial phase). Four of the nontrivial SPT phases, i.e., to,t x ,t y ,t z 
have been constructed in spin-| ladder systems [13]. However, due to a limited Hilbert 
space, the other three nontrivial phases cannot be realized with spin-]; ladders. Here 
we will realize all the eight SPT phases in spin-1 ladder models. 

The phase transition between different SPT phases is an important issue. It is 
known that in spin-1 chain there is a second order phase transition between the Hal¬ 
dane phase and the trivial phase if the Hamiltonian has a continuous symmetry, e.g., 
U( 1) spin rotation symmetry[30]. On the other hand, if the symmetry group is dis¬ 
crete, only first order phase transitions between different SPT phases were observed in 
the literature [1,3, 11]. In the following we study the transitions between different SPT 
phases and our results strongly support the following statement: in ID a direct contin- 
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uous transition between two SPT phases cannot take place if the symmetry group is 
discrete throughout. We will discuss about the special case where emergent continuous 
symmetry protects the continuous transition between two SPT phases (see [31]). 

A typical phase diagram of our model is shown in Fig. 1, where 0 1 <f> are the pa¬ 
rameters of the Hamiltonian [see Eq.(6)] and there are four SPT phases, one is trivial 
and the other three are nontrivial. Any two of these SPT phases are separated by at 
least one symmetry breaking phase. Direct phase transitions must be first order except 
that the symmetry group is enhanced to a continuous one. Our study will shed light on 
phase transitions between SPT phases with different symmetry groups and in higher 
dimensions. 

The remaining part of the paper is organized as follows. In Sec. 2, we introduce the 
linear and projective representations of the symmetry group D 2 x a, and give the active 
operators for each SPT phase. Further, we provide the Hamiltonian for each SPT phase. 
In Sec. 3, we study the phase transitions between different SPT phases, and some of 
the results are explained in Sec. 4 through nonlinear sigma model effective held theory. 
Sec. 5 is devoted to the conclusion. 

2. Parent Hamiltonian for each SPT phase 

Here we consider S = 1 spin ladders respecting D 2 x a symmetry. The discrete 
Abelian group D 2 x a has eight one dimensional linear representations (see Table A. 1 
of Appendix A). The two 5 = 1 spins on each rung span a Hilbert space 1 0 1 = 
0 © 1 © 2. The bases with total spin 0 and 2 are symmetric under P, they respect 
A g , Big, B 2g , B^g representations under group operation. The bases with total spin 
St = 1 are antisymmetric under P and respect Il\ „, B 2u , B 3u representations under 
the group operation . 

Z ?2 x (7 has eight classes of projective representations [11], the trivial projective rep¬ 
resentation is one dimensional, and the other seven nontrivial ones are two dimensional, 
as shown in Table A.2 of Appendix A. The projective representations correspond to 
the edge spins of the SPT phases. Due to their different edge states, these SPT phases 
can be completely distinguished by their different active operators, namely, by their 
different response to various external probing fields [11, 13]. 

In the following, we give the Hamiltonians of these SPT phases (the details of cal¬ 
culation are given in Appendix B). It has been shown in Ref. [13] that the t x ,t y ,t z 
phases can be obtained from the to phase via onsite unitary transformations (which do 
not commute with the symmetry group) since their Hamiltonians and active operators 
can be transformed into each other by these unitary transformations. However, the 
remaining three SPT phases t xy . t xz , t yz cannot be obtained by onsite unitary trans¬ 
formations from the Haldane phase to. In this sense, these three phases are more exotic 
comparing to the t x , t y . t z phases. If we collect the phases related by onsite unitary 
transformations into the same set, then the seven nontrivial SPT phases are separated 
into two disconnected sets, as shown in Fig.2. In the following discussion, we denote 
U\ (0, m ) as a rotation of angle 6 along ///-direction for the spins on one leg, and denote 
U(0, to) as a spin rotation on both legs. 

2.1. Hamiltonian for the trivial phase 

The trivial phase can be simply realized by two coupled antiferromagnetic Heisen¬ 
berg chains. The Hamiltonian is given as 

-^Trivial — J\ ^ ^ * S T) j-(-i © J 2 ^ B] A * S 2 ; i, (1) 

r,i i 
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Figure 2: Schematic figure for spin-1 ladder system and the relations between nontrivial SPT phases. In (b), 
the rotation only acts on one chain, while in (c) the U( 77 ) operation acts on both chains. 


where r = 1,2 labels the two legs, J\ > 0 is the intra-chain coupling and J 2 is the 
inter-chain coupling. 

Above model has SO( 3) x a symmetry. If J 2 > 0, the system has no edge states; 
when J 2 < 0, the system has spin-1 edge states. In both cases, the ground state is a 
trivial SPT state [32, 33], If the symmetry breaks down to D 2 x <7 by adding some 
anisotropic interactions, then generally there are no edge states no matter the inter¬ 
chain coupling is ferromagnetic or antiferromagnetic. 

2.2. Hamiltonians for the to,t x ,t y , t z phases 

The to phase is the analog of the Haldane phase of spin-1 antiferromagnetic chain 
since they both have 2-fold degenerate edge states and their edge states respond to 
external magnetic field in the same way. The Hamiltonian for to phase contains two 
parts, 

H 0 = Y^{H°^ 1+ H9), ( 2 ) 

i 

where Hf i+1 represents rung-rung coupling between rungs i,i + 1 and if' represents 
intra-rung coupling, 

=Y, Jo [ S -.‘ • S ^+! + A)(Sr,i • S T , I+ l) 2 ] 

r 

+ Jr [S T ,i ■ S w + 01 (Sr,i • Sr,i+i) 2 ] , (3) 

H? = Ji [Si,i-S 2l < + 02(S 1 , i -S2,i) 2 ] , 


4 


































with t = 1,2 and f = 3 — r label the two legs of the ladder. With parameters 
Jo = |,A) = 5 , J\ = |,/3i = J 2 = |,/32 = 1, above Hamiltonian H 0 is 
frustration free and its exactly solvable ground state has spin-1/2 edge states. The J\ 
term represents the diagonal interaction and can adiabatically be reduced to 0 without 
phase transition. Again, above model has an enlarged 50(3) x a symmetry which can 
be broken down to D 2 x a by introducing anisotropic interactions. 

The Hamiltonian H x {H y , H z ) of the t x (t y ,t z ) phase can be obtained from Ho by 
spin rotation of 7 r angle along x(y, z) direction for the spins on one of the chains. H x 
also contains two parts, 


H x = E (% 1 + ^), 


(4) 


H: 


i,i+1 


H? = 




Jo [S T ,i • S Tj i + i + Po(S T> i ■ S r j_|_i) 2 ] 


_L T lYO* Q x qV qV q z \ 

i D T,i D T,i+l °T,i D f,i+l) 

,/D I QX QX _ qV qV _ qZ qZ \2l 

_ rPlt* :, r,i 0 f,i+l J T,i°r,i+l J T,i D f ,i+l ) Ji 

T f/ q x qx qy qy qZ qZ \ 

I a f qx qx qy qy qZ qZ \2l 

x~P2\‘~ > l,i D 2,i J l,i J 2,i J) 


with (here m = x, y , z) being the spin -1 operators at the rth chain and ?'th rung. 
Similarly, we can obtain the Hamiltonians for t y and t z phases. 


2.3. Hamiltonians for the t xy ,t xz , t yz phases 

The Hamiltonians for t xy , t xz . t yz phases cannot be obtained by simple onsite uni¬ 
tary transformation of H (> . However, the Hamiltonians and the active operators of these 
three phases can be transformed into each other by onsite unitary transformations (see 
Fig. 2). For each phase, one exactly solvable Hamiltonian was constructed (see ap¬ 
pendix for details). For example, the Hamiltonian for the t xz phase is given by 

H xz = H 2 ^ + H 2 f 2 + Hf. +1 + Hl +1 (5) 

i 

1 r2 1 

where H- represents on-site term, H i ’ i+1 represents two-body rung-rung coupling and 
H 2 ' 2 represents intra-rung coupling, Hf i+1 and Hf i+1 represent three-body and four- 
body interaction. Due to its complexity, we leave the full expression of (5) to Appendix 
B. 

This Hamiltonian (5) has an enlarged symmetry. For instance, the Z 2 subgroup 
generated by R y becomes a Z 4 group with the generator U(^,y). (5) is also invariant 
under the continuous transformations (7(1) x (7(1), where the first (7(1) is generated 
by J2i( S i,i) 2 ~ ( S 2 ,i) 2 and the second (7(1) is generated by X7(5jF ) 2 - (S% ti ) 2 . 

The Hamiltonians of t xy and t yz phases can be obtained from on-site unitary trans¬ 
formations shown in Fig. 2 and will not be shown here. Contrary to the to (t x . t y . t z ) 
phase, the Hamiltonian for the t xz (t xy , t yz ) phase cannot be adiabatically simplified 
into two-body interactions. That is to say, the three-body and four-body interactions 
are important to realize the t xz (t xy , t yz ) SPT phase. 


3. Phase transitions between SPT phases 

In the previous section we give the Hamiltonians of the SPT phases. Here we study 
possible phase transitions between different SPT phases. To this end, we first consider 
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(d) 


Figure 3: 2D phase diagram with three SPT phases t xz ,to,t y . 6 = 0 is for H xz . 6 = 0 = 0(^)is 

for Ho(Hy). (a), entanglement spectrum gap A p = pi — p 2 , where pi and p 2 are two maximal Schmidt 
eigenvalues of the ladder; (b), entanglement entropy S = — trplnp, where p is reduced density matrix of 
the ladder; (c)(d), local magnetization order \{S^ ^)| and |(S'| t i )|. All the magnetization orders are in units 
of h, and the same for the following figures. 


the Hamiltonian 


H{9, </>) = cos 6H XZ + sm9(cos<j)H 0 + sin 4>H y ) (6) 

which connects t xz , to, t y phases through parameters 9, <f>. This model has three non¬ 
trivial SPT phases and the phase diagram can give a rough picture of phase transitions 
between SPT phases. The calculation is carried out using infinite time-evolving block 
decimation (iTEBD) method [34] with virtual dimension D = 40 (and the same be¬ 
low). 

The 2D phase diagram for Hamiltonian ( 6 ) is shown in Fig. 3. Three nontrivial SPT 
phases t xz ,to, t y manifest themselves in Fig.3(a), which are characterized by vanishing 
of entanglement spectrum gap and no magnetic orders. A trivial SPT phase locates near 
the center of the phase diagram, which has a large entanglement spectrum gap but no 
magnetic orders. The wave function for this trivial SPT phase can be adiabatically 
connected to a direct product state 

|trivial) y = 0 ^ (|1,0) - |0,1) + |0, —1) - | - 1,0 ))*, 

i 

which is the eigenstate of Sf + = + S^ i with the eigenvalue 0. Most of the 

phase transitions seem to be second order owing to the large entanglement entropy [see 
Fig.3(b)]. 
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(a) (b) 


Figure 4: Phase diagram with a direct first order phase transition between t xy and t X z phases. A = 1 
is the point for H xy , and A = 0 is for H xz . The entanglement spectrum gap A p, A p = 0 means the 
entanglement spectrum is doubly degenerate. The SPT phases are characterized by double degeneracy of 
entanglement spectrum and vanishing of magnetic order. 


In above phase diagram (also see Fig. 1), the SPT phases are separated by several 
symmetry breaking phases and never touch each other. That is to say, direct phase 
transitions DO NOT take place between any two SPT phases. The symmetry breaking 
phases have either ( S v± ) or (S z± ) order, as shown in Fig. 3(c) and (d) [35]. 

To see generally if there exist direct phase transitions between any two of the SPT 
phases, we connect the Hamiltonians H a ,H b of two SPT phases t a , t b via a single 
parameter, 

H(X) = \H a + (1 - \)H b 

and study the corresponding phase diagram. It turns out that two SPT phases may 
be separated by: (A) a direct first order phase transition; (B) a trivial phase; (C) a 
symmetry breaking phase; (D) a trivial phase and symmetry breaking phases. We will 
give an example for each case. 

3.1. Direct first order phase transitions 

The t xy phase can be changed into t xz phase via a first order phase transition. To 
see this, we interpolate a Hamiltonian path between H xy and H xz , 

H = A H xy + (1 - A )H XZ 

and study its phase diagram. This model has an enlarged global symmetry (7(1) x 
(7(1), where the first (7(1) is generated by fij 2 — (Sf J 2 and the second (7(1) is 

generated by £7(5 ^) 2 - (S 1 ^) 2 . 

The results are shown in Fig.4. The ground state energy E{ A) and its derivatives are 
shown in Fig.4(a). There is a singularity in first order derivative of the energy. Further¬ 
more, no magnetic orders are found and the entanglement spectrum is always two-fold 
degenerate except the transition point [see Fig.4(b)]. Above information shows that 
there are only two SPT phases t xy and t xz , and the phase transition between them is of 
first order. 

Direct first order phase transitions can also be found between t xy and t yz phases, 
or between t xz and t yz phases, or between any two of t x ,t y , t z phases. 
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(a) (b) 


Figure 5: Phase transitions interpolating to and t x . (a), energy and its derivatives; (b), the local order is 
|(5^ + )| = (—1 ) l (Sf i + 5*2 *)• Between two SPT phases there is a trivial phase characterized by a large 
entanglement spectrum gap A p and vanishing magnetic orders. 


3.2. Phase transition through a trivial phase 

Between to and t x phases we did not find a direct phase transition in the phase 
diagram of the model 

H = XH 0 + (1 — A )H X . 

This model has an enlarged (accidental) global 17(1) symmetry generated by S“ + . 

It turns out that the two phases to and t x are separated by a trivial SPT phase, as 
shown in Fig. 5. The trivial phase can be adiabatically connected to the direct product 
state 

|trivial) x = 0 t (|l, 0) - |0, 1 ) - |0, -1) + | - 1,0)) i . 

i 

At each rung above state is the eigenstate of Sf + = S'f + Sfy t with the eigenvalue 0 
[36], 

The phase diagram in Fig. 5 shows that both to and t x phases are narrow regions 
sandwiched by a trivial phase and a stripe Neel ordered phase with nonvanishing 
\(Sp}\. If the accidental 17(1) symmetry is removed by some perturbations, then 
an extra symmetry breaking phase will appear to separate the trivial phase and each 
SPT phase(see Sec. 3. 3). The physical reason behind this will be discussed in Sec. 4. 

3.3. Phase transition through a symmetry breaking phase 
In studying the model 


H = XH 0 + (1 - X)H xz , 

we find different symmetry breaking phases separating the trivial SPT phase and t 0 {t xz ) 
phases as shown in Fig. 6 . 

In Fig. 6 , there is a small S z+ order coexisting with the S x+ order at the same 
parameter region with | (S? + ) | = (—l)* + 1 (5f, + Sf This order is relatively small 
with the largest value 1.4 x 10 -2 , so we did not show it in the figure. 

The trivial phase in Fig .6 can be adiabatically connected to the direct product state 

| trivial) y = (g) \ (|1,0) - |0,1) + |0, -1) - | - 1,0)),. 
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(a) 


(b) 


Figure 6: Phase diagram for to and t xz phase transition through a trivial phase and two symmetry breaking 
phases. Between the trivial SPT phase and t X z phase, there is an antiferromagnetic ordered phase where the 
order is |(S'f - )| = (—1 ) l (Sf i — S% ^). While between the trivial SPT phase and to phase, an antiferro¬ 
magnetic ordered phase exists in which the order is 1(5'®"*’) | = (—l) 4 (Sf i + S% *)• The region A < —0.1 
is a new symmetry breaking phase with order | | = (— l) 1 {Sf\ + S^j)■ 


At each rung above state is the eigenstate of Sf + = Sf i + S% i with the eigenvalue 0. 

The phase diagram in Fig. 6 tells us two pieces of information: (1) the non-trivial 
SPT phases are sandwiched by symmetry breaking phases; (2) two SPT phases (in¬ 
cluding trivial SPT phase) are separated by a symmetry breaking phase. 

Now we study another example, 

H = \H X + (1 — A )H XZ 
which connects the t x and t xz phases. 

The phase diagram is very similar to the previous one. As shown in Fig. 7, the 
intermediate trivial phase can be adiabatically connected to the direct product state 

|trivial)= 0 \ (| 1 , 0 ) + | 0 , 1 ) - |0, —1) - | - 1 , 0 )),;. 


At each rung above state is the eigenstate of Sf = Sf i — t with the eigenvalue 
0. The sharp peak of order parameter ((S’? - )! (|(S'f z,+ )|) indicates that there is a 
symmetry breaking phase separating the trivial SPT phase and t x (t xz ) phase. 

Once again, the phase diagram in Fig. 7 shows that generally two SPT phases are 
separated by a symmetry breaking phase. Similar results can also be obtained between 
t x and t y phases, or between t y and t xz phases. 

3.4. Summary of phase transitions between SPT phases 

According to the numerical results we obtained, two different SPT phases are sepa¬ 
rated either by a first order phase transition or by intermediate trivial or/and symmetry 
breaking phases. Direct second order phase transitions between nontrivial SPT phases 
are never observed. 

Above results are reasonable for discrete symmetry group protected topological 
phases. Notice that nontrivial SPT phases are disordered by strong quantum fluctua¬ 
tions (also see Sec. 4), and the quantum fluctuations are resultant from competition 
between different classical orders. Since discrete symmetry can be spontaneously bro¬ 
ken in ID, SPT phases are always sandwiched by symmetry breaking phases if the 
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Figure 7: Phase diagram for t x and t xz phase transition through two symmetry breaking phases and a trivial 
phase. Between the two nontrivial SPT phases, there is a trivial phase and two antiferromagnetic ordered 



(S? )| order. When A < —0.15, the antiferromagnetic ordered phase with the magnetic order \(S y )| 
appears. 

symmetry is not enhanced [13]. This explains why nontrivial SPT phases are always 
separated by one or several ordered phase. A direct first order transition between two 
SPT phases is possible but is not interesting. Our results suggest that in ID generally 
different nontrivial SPT phases protected by discrete symmetry DO NOT have direct 
continuous phase transitions if the symmetry is not enhanced to a continuous one. This 
is an important message of the present paper. 

However, if the symmetry group which is broken in the intermediate ordered phase 
is enhanced to a continuous one, then the intermediate ordered phase will vanish in 
ID owning to Mermin-Wagner theorem [37]. In that case, the two SPT phases have 
direct second order transition. That is to say, the continuous phase transition between 
different SPT phases (including the trivial SPT phase) is ‘protected’ by continuous 
symmetry. This phenomenon has been observed in previous numerical studies [1,3, 
13,38]. 

4. Effective field theory of SPT phases 

In the following we will interpret some results of previous sections from low energy 
effective field theory. We will start from a single spin-1 chain model, discuss possible 
phases and phase transitions, then we couple two chains to form ladder models [39]. 

4.1. Nonlinear Sigma Model description of a single chain 

SO( 3) symmetry. For a single spin-1 chain Heisenberg model H = S, • 
Sj+j, it was shown that the low energy effective field theory is the topological SO (3) 
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nonlinear sigma model (NLSM) [2] 



( 7 ) 


where 9 = 2 ttS with S = 1. Denoting |n) as the spin-1 coherent state, then we can 
define Berry connection A t , = (n|// ; , |n) and the Berry curvature F = d T A x — d x A T . 
The //-term can also be written as S$ = i ® ff drdxF. 

The dynamic term ,SVi is unimportant since under renormalization group (RG) the 
coupling constant g flows to infinity. The //-term is very important, from which we can 
read off the ground state wave function under periodic boundary condition (here we 
ignore the dynamic term). 



( 8 ) 


X 


Under open boundary condition (supposing the temporal boundary condition is peri¬ 
odic), the //-term 



effectively behaves like Berry phase terms of spin-1/2 particle living on each boundary, 
which explains the existence of spin-1/2 edge states of the Haldane phase [40], 

0(2) symmetry. Now if we add an anisotropy term D(Sf) 2 to the Hamilto¬ 
nian, then the symmetry of the chain becomes 0(2) = (7(1) x Z 2 , where (7(1) is 
continuous spin rotation along 7-direction and Z> is generated by a spin rotation of 7r 
along ./'-direction. In the low energy effective theory, the dynamic term becomes 



In the large D —> +oo limit, the vector n lies in the //'//-plain and the integrand of the 
//-term vanishes. The resultant ground state wave function is a trivial product state 



where |0) is the 0-component of spin-1. 

With the decreasing of the value of D, the ^-component of n increases. After a 
critical point D c the //-term plays an important role and the system enters the Haldane 
phase. If we ignore the dynamic term, then the ground state is given by (8) under 
periodic boundary condition. At the transition point D c , the fluctuation in the 9 -term 
is so strong such that the bulk gap closes. 

When D further decreases, z becomes easy axes as D < 0 and at some critical 
D' c the systems form antiferromagnetic order (S z ) \, which spontaneously breaks the 


11 


Z -2 symmetry. So the Haldane phase is sandwiched by a trivial phase and a symmetry 
breaking phase. 

Z 2 x Z -2 symmetry. Now we further introduce anisotropy into the Heisenberg 
interaction, then the dynamic term becomes 

S d = j j drdx^[J x (d^ x ) 2 + {d y n y ) 2 + (d y n z ) 2 } + Dn 2 z , 

and the symmetry of the chain reduces into Z 2 x Z>. We will see that the Haldane 
phase is wrapped by three symmetry breaking phases. 

In the large J x limit, the systems falls in a symmetry breaking phase M x and has an 
antiferromagnetic order | (S x )\. In the large negative D limit, the system is in M z phase 
which contains antiferromagnetic order |(S' Z )|. On the other hand, if D is positive but 
not too large (otherwise the system enters the trivial phase), then the vector n favors 
xy- plane, at the same time if ./,. < 1, then the vector n favors yz- plane. Resultantly 
within a certain region of D and J x , the system will fall in the M y phase and has 
antiferromagnetic order |(iS'j,)|. 

The Haldane phase locates near the center of the three ordered phases, where all the 
three order parameters are fluctuating such that n can be pointing to arbitrary directions. 
Thus the 0-term is non-vanishing and gives rise to a SPT phase. In other words, the 
Haldane phase is a consequence of the competition between the three orders. 

Especially, the previous direct transition from the large D trivial phase to the Hal¬ 
dane phase is split into two phase transitions, since there is an intermediate symmetry 
breaking phase M y (if J x < 1) or M x (if J x > 1). Owing to the decreasing of the 
symmetry, the anisotropy terms D and J x suppress the topological 0-term and trigger 
the symmetry breaking, which prevents the touching between the trivial phase and the 
Haldane phase except at the high symmetry point J x = 1. 


4.2. Nonlinear Sigma Model description for ladder SPT phases 

to phase. Now we consider two spin-1 Heisenberg chains coupled with inter-chain 
interaction Hf = J [4 - (S M x S 2ji ) 2 ] = J[(S M • S 2ji ) + (S M • S 2ji ) 2 ]. Now the 
system has <S'0(3) symmetry and the effective field theory is given by 


Z 


J Dti(x 1 t) exp(— S d — 


Se), 


S d 


S e 


j j dTdx^[(d y ni) 2 + (d y n 2 ) 2 } - J(n d x n 2 ) 2 , 
id J j drdx^-lni ■ {d T ri\ x d x ri \) 


+n 2 • (d T n 2 x d x n 2 )\. 


(9) 


where 9 = 2ir. 

In the limit J —> oo, the interchain term behaves as a projector such that at each 
rung ni and n 2 are orthogonal to each other, namely ni(x, r) • n 2 (x , r) = 0. Under 
this condition, the topological term is equal to 


S e = id' 


drdt 


—to • (d T m x d x m) 
47r 


where m — (nr + n 2 ) is a unit vector. Notice that 9' = s/29 = 2\/2 -k is not 
quantized. However, since 2n < 9' < 37 t, under RG 9' will flow to its fixed point 
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27r [41, 42, 43], As a result, the strongly coupled spin-1 ladder is similar to a single 
Haldane chain and has spin-1/2 edge states (i.e. the to phase). 

Trivial phase. Now suppose J is small such that ni and «2 are almost independent. 
Then we can treat the interchain coupling J term as a perturbation. Owing to the 9- 
term for both n t and n 2 , under open boundary condition, each chain has a spin- 1/2 
Berry phase term at the end. The interchain coupling —J(n\ x n 2 ) 2 ~ ■) (n,\ ■ n 2 ) 2 
will cause antiferromagnetic correlation between the edge states of the two chains, as 
a result the Berry phase terms tend to cancel each other such that effectively there is no 
net Berry phase at the boundarytwhich is equivalent to 9 = 0). In other words, the two 
spin- 1/2 edge spins of the two chains form a singlet and finally the system has no edge 
state. So the system is in a trivial phase. 

When increasing J, effectively the 9 varies from 0 to 27 t (here the 9 is the fixed 
point value after RG flow), so there will be a continuous phase transition (at 9 = ir) 
from the trivial phase to the Haldane phase. 

t x phase. The t x phase differs from the to phase by the interchain coupling. Since 
the Hamiltonians of the two phases are related by an onsite transformation Ui(tt,x), 
the effective field theory for the t x phase can be obtained from (9) by replacing n i with 
hi = (nf, — n\ , — raf) T , namely, replacing —J(n\ x n 2 ) 2 with x n 2 ) 2 and 

keeping the rest part unchanged. 

Similar to the previous discussion, it can be shown that the t x phase with large J 
also has two-fold degenerate edge states. On the other hand, small J results in a trivial 
phase (if the system has no translational symmetry, this trivial phase is the same phase 
with the preceding one). Increasing J to a critical value, the system will undergo a 
continuous phase transition from the trivial phase to the t x phase. 

Phase transitions from to to t x . Above we have illustrated a possible path from 
the to phase to the t x phase: 

to —t trivial —> t x 

by tuning J to zero then turning on J. Similar result can be obtained if the interchain 
coupling is parameterized by A 

— [AJ(ni x n 2 ) 2 + (1 - A)J(ni x n 2 ) 2 ]. 

A = 0 gives the t x phase and A = 1 is the t 0 phase. In the intermediate region near 
A = 1/2, the interchain coupling behaves like a projection operator locking n t and n 2 
into a scalar, so the resultant phase is trivial. 

Until now the system has (7(1) symmetry, so direct phase transitions take place 
between the nontrivial SPT phases and the trivial SPT phase. Even though, differ¬ 
ent nontrivial SPT phases are untouched. If we lower the symmetry by introducing 
anisotropic Heisenberg interaction, then similar to the discussion for a single chain, 
even the direct phase transition between nontrivial SPT phase and the trivial phase will 
be prohibited, and intermediate symmetry breaking phases will occur. 

In principle, the effective theory for the t xz ,t xy ,t yz phases can be obtained from 
their microscopic Hamiltonians. However, since the interactions contain quadrupole 
terms and many body interactions, the effective theories are not the same as usual 
topological NLSMs. We will not discuss about the details here. 

In above discussion, the symmetry operations acting on the fields n \. n 2 are the 
same for all SPT phases, as defined by the symmetry group. On the other hand, if 
we allow the symmetry operations to act differently for different SPT phases, then the 
effective field theory for all the SPT phases can be written in the same form, i.e. the 
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SO( 3) topological NLSM [9]. If we discretize space and time, the effective field theory 
of SPT phases can also be described by topological NLSMs with group cocycles [8]. 

5. Conclusion 

In summary, we have studied SPT phases protected by symmetry group D 2 x a. We 
realized all the SPT phases on spin-1 ladder models and studied the phase transitions 
between different SPT phases with iTEBD method. We did not find continuous phase 
transitions between any two SPT phases. Instead, two SPT phases are separated by 
either a first order transition or by one (or several) intermediate phase. We interpreted 
part of our results via topological nonlinear sigma model effective field theory. We 
further conjecture that in ID generally continuous phase transition cannot take place 
between SPT phases protected by discrete symmetry group, unless the symmetry is en¬ 
hanced to a continuous one, or a continuous symmetry emerges near the critical point in 
the low energy and long wavelength limit. If the critical point between two SPT phases 
exists, then the local order parameters corresponding to discrete symmetry breaking 
have power-law correlations, these quasi-long-range orders can be easily turned into 
long-range orders (unless the order is associated with continuous symmetry breaking) 
by perturbations preserving the discrete symmetry. Our results provide some hints for 
studying the transitions between SPT phases in higher dimensions [31, 44, 45]. 

We do not rule out the possibility that the Hamiltonian only has discrete symme¬ 
tries but at very low energy and small momentum a continuous symmetry emerges. 
That is to say, in some special fine tuned models, direct second order phase transitions 
may exist between SPT phases protected by discrete symmetries owning to the emer¬ 
gent continuous symmetry at the critical point [31]. However, if symmetry preserving 
perturbations are added to the Hamiltonian, the direct continuous transition should be 
unstable and may split into several transitions between which discrete symmetry break¬ 
ing phase appears. 

6. Acknowledgement 

We thank Li You, Hong-Hao Tu, R. Orus, Peng Ye, Zheng-Cheng Gu and Chao- 
Ming Jian for helpful discussions. JYC is supported by the MOST 2013CB922004 of 
the National Key Basic Research Program of China, and by NSFC (Nos. 91121005 and 
11374176). ZXL is supported by NSFC 11204149 and Tsinghua University Initiative 
Scientific Research Program. We acknowledge the computation time on the clusters in 
IASTU. 


Appendix A. linear and projective representations of the symmetry group 

In this appendix, we give the linear and projective representations of the symmetry 
group D 2 x a, as shown in Tab. A.l and Tab. A.2, respectively. 


Appendix B. constructing the hamiltonian for SPT phase 

In this appendix, we show the method by which we obtain the Hamiltonian for each 
SPT phase. The method is closely related to the way for obtaining Affleck-Kennedy- 
Lieb-Tasaki (AKLT) Hamiltonian [46]. Here we follow the method used in Ref. [13]. 
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The procedures contain three steps:(i) construct a matrix product state (MPS) wave 
function with given edge states which are described by the projective representations; 
(ii) construct the parent Hamiltonian for the MPS using projection operators; (iii) sim¬ 
plify the parent Hamiltonian by adiabatic deformations. The linear representations 
have been shown in Table A.l, and the projective representations in Table A.2. The 
Hilbert space of the direct product of two projective representations can be reduced to 
a direct sum of linear representations, which has already been listed in Ref.[ll], In 
the following, we will illustrate the method using the Hamiltonian of to phase as an 
example. 

For to phase, the edge states can be described by E 5 and E e representations. (E e 
is not shown explicitly in Table A.2. Both E e and E s can be found in Ref.[l 1].) In an 
ideal MPS, every rung is represented by a direct product of two projective representa¬ 
tions, which can be reduced to four linear representations. Here we choose E 5 0 Eq 
with the Clebsch-Gordan (CG) coefficients: 

E§ 0 Eq = A u 0 B\ u 0 B 2u 0 Bou 

C A “ = a x , C B = ia y , C B = <r s , = /; ' ^ ' ' 

where the CG coefficients have been chosen real. From the CG coefficients, we can 
write an ideal MPS which is invariant (up to a phase) under the symmetry group D 2 x er: 

\t) = ^ Tr{A m K..A mN )\m 1 ...m N ) (B.2) 

with A m = e l0m BC" 1 . Here B is the CG coefficient of decomposing the product 
representation £5 0 Eq into one-dimensional representation which was chosen to be 
C Au and e“ m can be absorbed into the spin bases. Thus we have the four bases on 
each rung: 


l</>i) = -^(l 1 ,- 1 )-|-M»; 1^2) = 11) 0) — 10,1); 
l</>3> = 10, -1) + | - 1,0 ); 1^4) = ^(|1, -1) - I - 1,1)). 


(B.3) 


Table A. 1: Linear representation of D 2 X a, P is the interchain permutation operator. Notice the notation 
= 0 1 d= O 2 is a single-rung operator, where 0\ is the operator for the first chain and O 2 is that for the 
second chain. 



R z P 

bases 

operators 

A 

1 1 1 

| 0 , 0 >,| 2 , 0 >,^(| 2 , 2 ) + | 2 ,- 2 » 

02+ q‘ 2 + 02 + 

°x ’ °y ’ °z 

Big 

1 -1 1 

i(| 2 , 2 )-| 2 ,- 2 » 

St 

B 2 g 

-1 -1 1 

i(| 2 ,l)-| 2 ,-l)) 

s+ 

Big 

-1 1 1 

i(| 2 ,l) + | 2 ,-l)) 

St 

■Ay, 

1 1 -1 


S't, Sf, S‘t 

B\u 

1 -1 -1 

|i,o) 

Sz 

b 2u 

-1 -1 -1 

^(| 1 , 1 ) + | 1 ,- 1 )) 

s t 

Bou 

-1 1 -1 


S J 
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Table A.2: Projective representations and the corresponding SPT phases of D 2 X a. The effective operators 
for the active operators are ((T x , cr y , cr z ). The active operators can split the degeneracy of ground states on 
open boundaries. 



R z 

R x 

p 

Active operators 

SPT phases 

Eg 

1 

1 

1 


I (trivial) 

Ei 

I 

icy z 

a y 

(St v . z ,S 7 ,St) 

txy 

E3 

cr z 

I 

ia y 

(Sxv, Z ’Sx’S+) 

ty Z 

e 5 

ia z 

& X 

I 

(S+, s+, St) 

to 

e 7 

or 

io z 

i<j x 

(S'i~, z ,S-,S+) 

txz 

Eg 

icr z 

OX 

i(j x 

(St, S-, S-) 

tx 

En 

ia z 

io x 

cy z 

(S-, S-, S+) 

t z 

E13 

ia z 

i&x 

iCTy 

(S-, St, S-) 

ty 



(a) (b) 


Figure B.8: Deform the exact Hamiltonian of to phase continuously with a parameter d. When d = 1, it 
is the exact Hamiltonian. When d = 0, we have dropped the three-body and four-body interaction terms, 
(a) Energy, first energy derivative and second energy derivative, (b) Largest of Schmidt eigenvalue and 
entanglement spectrum gap. The numerical method is iTEBD with virtual dimension D = 40. 


The orthonormal bases for the MPS state of two neighboring rungs i,i + 1 are: 


IV’l) — ® I0l)i+1 + I $ 3)1 ® |02)i+l); 

1 ^ 2 ) = 2 (l^l)i ® I h)i+l - 103 )i <8> |0l)i+l) 

I03> = \{\^>2)i ® |0l)i+l - |0l)i ® \<t>2)i+l) 

IT 4 1 

IlM = \ ® ^3)1+1 + -=\<t>l)i ® \<i>l)i+l - 7 # 3 )i ® |</> 2 )i+l). 

V 24 5 5 


(B.4) 


The parent Hamiltonian iJ 0 ,es = Ei H i = ~ E, Ea=i iV’aKV’a] has many terms. 
Hi = H 2 + d x (//•{ + H 4 ), where d = 1, //.j and H,\ represent three-body and 
four-body interaction, respectively. When we tune the parameter d from 1 to 0, there 
is no phase transition, shown in Fig. B.8. Further dropping the constant leads to the 
Hamiltonian in Eq.(2). 

The same method can be used for t xz phase, where we choose projective represen- 
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tation £7 < 8 > Eg as our starting point. 


£7 g) Eg, — Ag ® B 2g © A u ® B 2u , 

C A ° = a z ,C B2a = ia y ,C A “ = I, C B = a,; 


(B.5) 


and we choose 10, 0) as the representation basis for A g . Following the same procedure 
and dropping the constant, it leads to the Hamiltonian (5) with the expression: 


h xz = Y J H l+ #5+i + H ?' 2 + H li+ 1 + H t+i 

i 


(B. 6 ) 


r = £> [S v T ,S y T , +1 + (S?,<S?,<+1 + 5 'iS'i+i) 2 ] 

T 

+ VF 2 (5^5^ i+1 ) 2 + W3[(S^ ii+1 ) 2 + (S^+i) 2 ] 

+ W 4(^^ ii+1 + 5*^5 i+1 ) 2 + WsKS^+r) 2 + (5^5? i+1 ) 2 ], 


H 2 ’ 2 = W 6 (S M • S 2il ) 2 + WrKS ?,^,*) 2 + (S z hi S^) 2 - (Sf^.) 2 ], 


Hf ti+1 = C/r(S M • S 2jl ) 2 [(^ +1 ) 2 + (5" i+1 ) 2 ] + £1 [(Sy 2 + (5'| i j) 2 ] (S M+1 • S 2 , i+1 ) 2 
+ U 2 (S? ti S & 2 [(Si , +1 ) 2 + (S 2 * >1+1 ) 2 ] + U^Sl.S ^) 2 [(Sf ii+1 ) 2 + (Sf, i+1 ) 2 ] 

+ u 2 [(Si ,) 2 + (S^f] (Sl i+ 1 S z w f + U 2 [(Sf,*) 2 + (Si,) 2 ] (Si, + 1 Si , +1 ) 2 

\TT \(Q x Q x y \ Q z QV Z Q x v Q x qyz qz \fqy qy \] 

' u 3LV°l,i°2,2 ' °l,i D 2,i °l,i D 2,i)\ lD l,i+l °2,i+l)\ 

_ TT r ( QX qz J_ qz qX , qXy qyz , qyz qXy \ / q X Z , qXZ \1 

u ^l\°l i i D 2 t i ' 0 l,i 0 2,i ' °l,i°2,i) V°l,i+1 ■ °2,i+l)\ 

j_tt \(QV Qy \(q x Q x y q z qy z q x y qx qyz qz \i 

-r ^3LWl,i °2,i/l 0 l,*+l 0 2,t+l ' °l,i+l°2,i+l °l,i+l°2,i+l 0 l,i+l°2,i+l/J 

_ TT r ( q XZ 1 qXZ\(qX qZ . qZ qX , qXy qyz , qyz qXy \1 

^3Lv°l,i 1 0 2 ,iA°l,i+l° 2 ,i+l ■ °l,i+l 0 2 ,i+l ■ °l,i+l°2,i+l ' °l,i+l 0 2,i+l/J 

+ u, [susi ,) 2 - (Sl,) 2 si,] (s \, +1 - si ii+1 ) 

+ C / 4 [Sf.KSZ ,*) 2 + (S?^) 2 ^] (5f,l + i + Sfji+t) 

+ £ 4 ( 5 ^ - S^) [S?, j+ i(SZ, j+ i) 2 - (5f ii+1 ) 2 5| ii+1 ] 

+ t/ 4 (^ + [S% + 1 (S £ i+1 ) 2 + (Sl i+ 1 ) 2 S% +1 ] 
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Hii +1 = Wm ' S 2 ,i) 2 [(S x hi+1 Sli +1 ) 2 + (S?, j+1 S|, i+1 ) 2 - (S? ;i+1 S| ii+1 ) 2 ] 
+ M^,A*i) 2 + - (S? >i ^,i) 2 ](Si, j+ i • S 2 , i+1 ) 2 

i t/ r ( Q z Q x Q x Q z \2 I f qx qz qz qx \2] 

' v 2 l{ 0 l,i 0 2,i°l,i+l 0 2,i+lJ ' Wl,i 0 2,i°l,i+l 0 2,i+lJ J 

+ ^3(Si, I -S 2>! ) 2 (S lii+1 -S 2 , i+1 ) 2 

I T / f QV QV I QXZQXZ\fQV ny _|_ QXZ QXZ \ 

' ' 0 l,i°2,i) V°l,i+l°2,i+l ' °l,i+l°2,i+l ) 

I T r ( QV qxz qxz ay \f Q y qxz qxz qy \ 

■T" v 4\°l,i°2,i °l,i D 2,i)\ D l,i+l°2,i+l 0 l,i+l 0 2,i+l/ 

+ J2ms~,s^ + s:”s^) + + s%s%) 

T 

+ VriS^iS^) 2 + (S^) 2 S^ + y 8 (S-(%) 2 + (S^) 2 S^)Y 

f qx qz . qyz qXy \ 

Wr,i+l°f,i+l "r °r,i+l°r,i+iy 

+ [V 5 (S^ - s^sf,) + v 6 (S^s;y - s%s* fii ) 

+ VrdS^) 2 ^ - Sf, { Sf,)f + VsdS^fSl, - S^S^) 2 )}- 

/ qx qXy _ qyz qz \ 

Wr,i+1 ^Tji+l ^r,i+l^r,i+l/ 

I [T/ / q x y qx qz qyz\ _i T/ / qyz qz qx qXy \ 

+ _ + ' 10 ~ 

+ Vnis^is^) 2 (S^) 2 sy + vMsf/s^) 2 - (S^fsy]- 

(S y Tti+1 (Sl i+1 ) 2 - (S z Tti+1 ) 2 S y w ) 

+ [V 9 (s^s z fti + s%s%) + v 10 (S;,s*, + s^sgi) 

+ VMS^Sl,) 2 + + V 12 (S%(S^) 2 + (S^fs?*)}- 

(( q x \2 qxz J_ qxz ( QZ \2\ 

U°r,2+1/ °f,t+l T ) 


where the parameters are given as: D = 4^, W\ = W 4 = — W 2 = —J, W 3 = 
f s ,W 5 = W 6 = -if ,W 7 = lU, = \,U 2 = | ,U 3 = y 2 .,U 4 = |, V 1 = 
l,V 2 = V3 = -if,V 4 =V 6 = ±,V 5 = ±,V r = V 10 = ±,V S = V 9 = 
^g,V\i = jg,V ± 2 = J 2 , and r = 1,2 ,r = 3 — r label the two chains, S mn = 
S m S n + S n S m (m, n = x,y, z, with m ^ re). 

The active operators have been discussed in Ref. [11] and will not be repeated here. 
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